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ABSTRACT 
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CONFIGURATIONS UNDER NATURAL CIRCULATION 
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Analytic comparisons of the performance of LMFBR open- 
hexagonal, triangular, and square lattice configurations 
under naturai circulation flow conditions were made. 
Lattices having the same volumetric compositions were 
compared with respect to reactivity change upon sodium 
voiding, pressure drop, and linear power output. 


A brief study of the relative advantage of tne lattices 
from the svandpoint of reactivity change upon sodium voiding 
of a coolant channel showed that the open-hexagonal lattice 
was preferable because of its larger coolant-channel radius. 


As regards hydraulic and thermal aspects, it was found 
that core pressure drops and linear power outputs were 
eameavelcate tor the three lattices under botn natural 
circulation and forced convection, for both laminar and 
turbulent flow regimes. 


Consideration was given to developing more accurate 
we waetona tor the hydraulic characteristics of the open- 
hexagonal lattice by modifying the traditional hydraulic 
diameter concept. Finally, some limitations of the open- 
iesaoonelelayvyblce as compared to the triangular lattice were: 
discussed and necessary areas for future work were identified. 
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Chapter 1 


TMeGoOdue.. on 

1.14 Background 

In liguid-metal cooled fast-breeder reactors an area 
meppervicular concern is that of loss of heat removal 
capability and subsequent coolant boiling. The designs 
for most large LMFBR's of the future include the capability 
for the removal of generated heat by natural circulation in 
the event of loss cf pumping power. Most of these designs 


here Dased on the use of triangular pitch fuel pin lattices. 


Hee Objectives 
Tne objective of this thesis is to determine the 

relative performance of different LMFBR lattice configura- 
vions under natural circulation with respect to the power 
output which could be attained during normal operation and 
With respect to safety considerations in the event of loss 
of primary pumping power in a pumped circulation system. 
Figures 1-1 and 1-2 describe the lattice configurations 


under consideration, 


ome Da scussion 
Advantages of natural circulation operation include: 
1, At full power operaticn primary pumps would 
be eliminated, a factor which would reduce cost 
of hardware and eliminate the concern of a loss 


of pumping power. 
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2. At full pover operation a natural circulation 
system would, to some extent, be self-regulating; 
for a given channel, as power increases/decreases, 
the core inlet and outlet temperature difference 
and the coolant mass flow rate would increase/ 
decrease. The parallel channel problem requires 
further investigation, however. 

3. In a pumped system natural circulation would 
serve as a safety measure in that the loss of 
primary pumping power would not result in coolant 
boiling and burnout since a sufficiently large 
natural flow could carry away the heat guickly 
enough. 

Disadvantages of natural circulation operation include: 
1. The plant will have to operate at a lover 
power output than would a plant with primary 
coolant pumps if design variables cannot 
adequately be optimized within various constraints. 
2. The design of the system would necessitate 
placing the intermediate heat exchanger far above 
the core in order to achieve the necessary heignt 


between the core and the heat exchanger for a 


desired thermal driving potential. In this 
position the intermediate heat exchanger could 
complicate access to the core from above and ang 
necessitate a larger Gealelcr) outer containment 


structure, 


eos 
1,4 Choice of Design Constraints and Variables 

The loop-type system was chosen as the primary coolant 
system. Since the reactor system needs a high thermal 
Griving head in order tc attain the desired coclant flow 
rate, the intermediate heat exchanger needs to be above the 
core at a height which cannot practically be designed into 
a pot-type system. However, the pot-type system, through 
@eewunermal inertia of its large pool of sodium, provides 
ample margin for core cooling following loss of pumping 
power, 

Sodium is by far the coolant of choice for recent LMFBR 
Gesi2ns. Thus, in this study liquid sodium was chosen as 
the coolant. 

Wie ctype of spacer system chosen for the open-hexagonal 
lattice was the grid system. Wire-wrap spacers would 
probably be more practical onemicaliy =), but with the 
center of the open-hexagonal lattice ere each rod would 
have support only from three wires in one complete spiral 
instead of support from six wires as in a triangular pitch. 
iiewerria spacer was chosen in order that the open~hexagonal | 
system would have the most and sufficient support against 
vibration, 

in the comparison of the three lattices several 
Variables will be examined. Under natural circulation 
the linear power output depends on a combination of core 


temperature rise, fuel pin diameter and pitch, fuel pin 





oon ee 
i neo, and the effective heat exchanger height to primary 
Coolant system length ratio. Combinations of the variables 
Mmemebe Ciscussed with the use of typical values of current 


LMFBR fluid and mechanical parameters. 


£25 Synopsis 

ive GOmparison of the LMPBR lattices will first be 
examined from the standpoint of the additional leakage 
meaeuiviuy resulting from the neterogeneity of the core 
using the different lattices. An examination of the 
relative mass flow rates, core pressure drops, and linear 
power outputs will then be made followed by an examole 
calculation of these quantities. In addition a more appli- 
cable expression for hydraulic diameter will be sought for 
an open-hexagonal lattice. Finally, conelusions and recom- 


mendations will be proposed. 
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Chapter 2 


Bee ncaculvity Eifect Components in a Sodium Voided Coolant 
Channel _ 


When a coolant channel is voided of sodium, a change in 
Peactivyity occurs. This reactivity change consists of three 
imeeoteoraer Components: a change in the number of neutrons 
captured, a shift in the neutron spectrum, and leakage of 
Memerons, Devaited discussion of the total sodium-void 
effect may be found in publications by Hummel and olsasel=! 
and Monroe!) , | 

Mose Cavensacvions of the leakage component of the change 
in reactivity are based on the assumption that the reactor is 
homogeneous. in this chapter a rough calculation of the 
additional leakage reactivity effect due to the heterogeneity 
meee iewcore will be made. For a more detailed discussion of 
the leakage component of sodium-void reactivity in hetero- 
Memeota cOress,) une reader is refered to a study made by 
Meeepelie’ whose method was based on the work of Beaostggh 
The leakage component of the overall reactivity change will 
be compared in the open-hexagonal and triangular lattices 
showing that the larger hole of the @pen=hexagonal cell 
results in a greater negative reactivity than does the hole 
Mmcweomurlenoular cell, The calculations used in the com- 
Dearison are based on formulations derived in Behren's 
original paper "The Effect of Holes in a Reacting Material 


w (6) 


on the Passage of Neulrons 





a SE 


2.2c Leakage Reactivity Calculations ae Open-hexagonal 
and Triangular Lattices 


Comparing Behren's Eq. (8) and Eq. (10) for closely 
spaced, interacting holes in a medium having a long mean 


free path we have 


2 
i Qr, 
=) = | + ode & + pases j (2-1) 
L y,\ 
O 
where 
) is the volume ratio of the holes to the Weteere kel: 


at che reactor, 


(8) is a function of the shape of the hole, 


r Poeoomiicane Gee pauhn in une solid, 

Phy = 2dy is the nuclear hydraulic radius of the eo 

y is the volume of the cell of material surrounding 
Seolnl Jeelior 

S is the surface area of each hole, 

Le is the mean square direct distance travelled by a 
neutron in the course of N collisions in a homo- 
geneous medium, 

ne is the mean square direct distance travelled by a 
neutron in the course of N collisions in a medium 
with holes, 


iierchiagdee in reactivity is given by 





Ok OD : 
k a ) (Kee ) oa ( 


The definition of diffusion length gives 
be pale 


_-<+ 





D 
= (LE) ele (2-3) 
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since 
O O O 
_ *aVeuel _ +a euel _ 2 
g, - ~Sfuel - —“afuel a(t) 
Veot Cia ncie) 62+ ¥) 


Equating Eqs. (2-1) and (2-3) 


(1 + i (1 + 9)? + hye 
Do d 
or 
Dea Dat tp) Phy PO , (255) 
e 3(1 +o) 


mae change in the diffusion coefficient is 


Qr, y 
6D =D - Ditom = __ hy" hole (226) 
om 
S 
where 
Diy is the homogeneous value of the diffusion 
om ae 
coefficient 
d 
Viole = aren is the volume fraction of holes 


For a first approximation we will assume that the cells 
im our theoretical core are cylindrical. The volume of each 
Sylandrical hole will be set Sate to the actual volume of 
Seebanu chenmel associated with a cell, and the volume of 
fuel material surrounding each hole will be set equal to the 
Towa levyolumenot fWel material associated with a cell, Refer 


Pome we kee for Fae Unit cell configurations chosen. Note 





-138- 
meat one cells are basec on the coolant channel and not on 
eeivel rod as is the more common oo any en, Vince tins 
approximation is made, the nuclear hydraulic radius, hy? 
is the same as the actual radius of the equivalent cylin- 
@erecal coolant cell. 

| Rotor eGo Mie, l.),a tor a definition of parameters 


needed to find the hydraulic radius of the open-hexagonal 


cell: 


oie = eae 9 4isnenme 


(21) 


Refer to Fig. 1.1.b for the hydraulic radius of the 


Welangular lattice: 


hole Ae ot 7 Maen 2 


2 ae 
rua )e Dy. as BP 5) 
ie 
A 2 2 hye 
(ie, = ley a ae (2-8) 
lapigue ay @ a G a P 
Since p, = (2/3)1/"p, then 
3 He 1 
2 2 B/N 
: (rip al Pr i 2 a Arias e (2-9) 





lS 
Inserting Eq. (2-6) into Eq. (2.2): 


ok - _(cuhy holey(** 7 ) 


= : (2-10) 
K 20 


k 3D om 


Using the following typical values for an open-hexagonal 


and a triangular lattice: 


Q = 4/3 (as tabulated in ref. (6)), 
Yrote = 1 = %ruer = 9-528 
Py =O ae LN, 
@, = cle ete 
d, = 0.300 in, 
= ) oo } Hf 
Dom + om imo ei eehOr- a. VOolGed core, 
ee = 1.5), 


ies ravio Of reactivity changes is 


(ok/k),, ; 
(Ok/k), 


solving for the reactivity changes: 


| 
i 


(ok/k), = -0,0056 = -$1.70 


(ok/k),, On Olle 


II 


-$3.40 


Pein p= 0.0033 = $1.00. 


mer a square lattice with 


Varna = (0.54 
OOO! Mi. 

qd, 3 2: 

p Om oo 11, 


S 





~OO- 


(ok/k) 
oS ieee 2 
( KM) og 

and. 


(Ok/k) = -$2,40 , 

From the above calculations we can see that the open- 
hexagonal lattice results in a higher negative reactivity 
from leakage in a voided coolant channel than do triangular 
or square lattices of the same volume fraction of fuel. 
Because of its higher negative reactivity the open-hexagonal 
lattice is preferable to both the triangular and the square 
lattices from che standpoint of keeping the effective 
imaerease an multiplication factor as low as possible upon 
feawum yOolcing of coolant channels. 

In the next chapter the relative fleco lon ray Cs, Core 
pressure drops, and linear power outputs of the three lat- 
tices will be compared. These comparisons will be used in 
Getermining which lattice is preferable from hydraulic and 


power removal aspects. 
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Chapter 3 


Expressions for and Comparison of Mass Flow Rates, Pressure 


| em a 


Drops and Linear Power Output Boe bistememe EMPBR Lattice 
Configurations 

3.1 Discussion 

As was seen in Chapter 2, the hexagonal lattice proved 
to ke more advantageous than a triangular or square lattice 
from a leakage reactivity standpoint. In this chapter the 
eee cmbattices Will >e compared with respect to some hydrau- 
lic ana power aspects to see if the larger hole of the open- 
hexagonal lattice renders it preferable to both the square 
and triangular lattices. In the present chapter the three 
Metices will be compared on an equivalent basis, namely, 
equal rod size and fuel volume fraction, with respect to 
pressure drops and linear power output. It will be assumed 
E@aeeali tluid properties are equal for the three lattices. 

As will be seen in the numerical examples of Chapter 4, 
the core pressure drops of the three lattice arrangements are 
essentially the same for equal fuel volume fraction and equal 
coolant velocities. The only differences in pressure drop 
are due to the differences in spacer arrangement of the 
different lattices, and these differences are small. MThere- 
fore, these small pressure drop differences can be neglected 
in a first order approximation of the relative linear power 


output for the three lattice arrangements. However, in the 


sODe 

calculation of the ASO lnee linear power output for any given 

mieerece, tne Spacer pressure Grop contribution to the total 

core pressure drop cannot be neglected since it constitutes 

as much as 45% of the total core pressure drop. 

Sze Procedure for Calculation of Core and Spacer Pressure 
Drop 


The total pressure loss, AP PaBough ene reactor Core, 


ie 
hence, through a unit cell flow path in the core, consists 
of three factors: pressure loss at bundle entrance and 

eat, AP.i frictional pressure loss without spacers, AP pi 


pressure loss due to spacers, AP Tt will be assumed that 


s° 
pressure loss at the bundle entrance and exit will be equal 
for all three lattices. This assumption is based on the fact 
that the entrance and exit openings are of equal total cross 
sectional area for the cores containing the different lattices 
and that the coolant velocities in the three cores are equal. 
Eyeusing equations Similar to Eqs. (3-1) and (3-2), one can 
find that the entrance and exit losses for the different cores 
are the same. 

Since there is no experimental data on pressure drops 
through open-hexagonal rod bundles with and without spacer 
devices known to the author, the approach for calculating 
the semi-empirical pressure drop through an open-hexagonal 


(7) 


lattice will be as follows. Tong's method for finding 
the theoretical pressure drop through a channel with Spacers 


is used to find theoretical values of AP. for the triangular 





aoe 
and open-hexagonal lattices. Then a AP. ioe OLUNG. fF Om 
experimental data gathered and correlated by de Stordeur ‘”? 
for triangular lattices. A ratio of experimental to theo- 
retical spacer pressure drops is ae obtained for trian- 
gular lattices. The theoretical pressure drop for the open- 
hexagonal lattice calculated from Tong's formulae are then 
multiplied by the ratio found above to obtain an appropriate 
semi-empirical spacer pressure loss for the open-hexagonal 
lattice. This approach thus uses the existing experimental 
data of cde Stordeur to normalize the calculated results. from 


Tong's method which is applicable, unlike de Stordeur's data, 


to both triangular and open-hexagonal lattices. 


3.3 Theoretical Core and Spacer Pressure Drop 

According to Tong "the pressure losses in a rod bundle 
fi weO the Lod Spacers are the form drag type pressure loss. 
This type of pressure drop can be calculated by using the 
pressure loss coefficients at sudden contraction, Kae and 
at sudden expansion, Kae as given by Kays and ation nu a el 
Values of K. and K, are Beopmocdnecd in Hig. 3.1. The pres- 
sure losses at sudden contractions and expansions are given 


as follows: 


For a sudden contraction 


2 2 
ep. 2 a 

AP, = ——(l-o') + ee per spacer. (s-a) 
256 ache 





“2 
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Fig. 3.1. Pressure loss coefficients at sudden 
COnbraccion ; K » and at sudden expansion, K as given 
by Kays and Londen 9 » Taken from Ref. (7), Fig, 1 
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For a sudden expansion 


VD eN5 
AP , =- —(l-o') - K.——  +=per spacer, (3-2) 
2g 29 


og = A,/A, is the ratio of the flow in the gridded region 
to the flow area away from the gridded region, 


is the velocity in the restricted cross-section 
region of the channel. 


Note that positive AP's are pressure losses. Refer to 
hit s.2 £LOr definition offlow areas. 


Therefore, 


AP hn = AP. = AP. ra 





1c (AP.,) - (8) 


where 


Momo micatLloOeGraticri nH ction factors in non- 
(10) 





Sime rrCcuUlamewEGecrlLecular Channels; see Fig. 3.3 
Substituting an expression for the friction pressure drop 
Eas aNet into Eq. (3-3) gives 
ree 
4olV, £ 


ae) + AP . + AP, : (3-4) 
De2zg fe 
e-"c 


AP a = AP. = fe 





Gixe 


Substituting expressions for AP, and AP, Peome bas.) (s-l)) and 


(3-2) into Eq. (3-4) 


4plV; £ p(A,/Ay) “Vi 
(——) + (K, + kK) 


De-Je eee 


(3-5) 
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Slo = 









— 


COoLane 


—————— A Spacer 


cupeanemmmammansnnseninly 


Region blocked 
by spacers 


iemeEeneouwmer Mow areas 14m Qa coolant cell. 


Ay is the area away from the spacer region, Ay is 
the free flow area in the spacer region, 


figs © 


uf OOO 


f/f cap 





Fig. 3.3. Fully developed turbulent flow parallel to 
Embaniewot scirculer tubes Of rods; points © based on friction 
eocificients caiculated by Deissler and Taylor‘ 11) at Re = 10°; 
Reynolds number influence is small, and Nusselt number 
behavior is virtually the same as friction behavior. Taken 
from Ref. (10), Fig. (9-10), which also shows that the ratio 
of Nusselt numbers is similar to the ratio of friction factors 


for non-metalic fluids. 


ae 
where 
Dy is the number of spacer grid sets in the rod bundle. 
Using 
fo = 0.046 Renee for turbulent flow, 


Eq. (3-5) becomes 
OP 2-0 Oo 1.3 





5 CORI Taney, f 
> 
A E De 9g fee 
e Cc cir 
p (A, /A,) Ve 
a (K + K Dr ; (3-6) 
2g e cg 
Cc 
where Va rom ounde: LOm.y —— w/PA, (ity ie ieeeand © 1S found 
ii Om 
2 g8o°Z,.D'**s"* 9g, 0.357 
Ww = a.) (ley ona) ier (A-8) 
20220 = Nc. 


The quantity w is the average mass flow rate for natural 
Circulation derived in Appendix A, ae is the number of effec- 


magesOds per cell and q = nail. 


3.4 Empirical Spacer Pressure Drop 


An empirical AP. for a triangular lattice will now be 


Bena using de Sterdeur's ‘®) meenod. According to de 
Stordeur 
C.VeSp 
1 
where 
Co Momigemspaces Grag coefficient 
V is the velocity in the spacer region (ft/sec) 


S 
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RA® 2 
S is the spacer projected frontal area (ft) 
g is the gravitational acceleration (£t/sec*) 


A is the flow area in the cell away from the 
gridded region (ft?) 


p is the density of the coolant (1b/£t*) 

AP. is the pressure loss due to one spacer (lb /£t7) 
Using appropriate values of Vor ono 5 Ay, 0, anda Co 
taken from Fig. 3.4, an empirical spacer pressure drop, AP as 
can be calculated from Eq. (3-7) for a triangular lattice with 
a honeycomb pin support grid. 

We can now obtain the ratio of the empirical spacer 
pressure drop to the theoretical spacer pressure aos how 
a triangular lattice with a honeycomb grid and set it equal 
menthe ratio of the empirical to theoretical spacer pressure 
drop for an open-hexagonal lattice: 


(AP. em] _ 


(AP) 


(AP) 
Sb em . (eae) 


sden}e L4?s) tntn 


An appropriate (AP.) on Gameoestounda from Eq. (3-8) for an 


open-hexagonal lattice after a (AP for the open-~hexagonal 


s) th 
lattice has been calculated using Tong's method. Equation 


(3-6) then becomes 


O20... 3 





weoo2 1p Vy % 
(APD ~ APE) * See en 
e Ie el ts 
(5) 
Cee) A(A, /A,) V 
S‘ em Lee 1 so 
(—_-—_———_) ecg 1 st Ko) Mg ° (3 Dy) 


(APS) en 2 2h. 
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Fig. 3.4. Drag coefficients for honeycomb grid spacer. 


(Dashed line is an extrapolated portion of the curve.) The 


quantity d., is the spacer width, °%) 


Sil 


3.5 Mass Flow Rate and Linear Power in Steady-State Thermal 
Ciitculatdaon 


In steady-state thermal circulation the thermal buoy- 


ancy driving force, AP equals the friction pressure drop, 


B? 


AP EOE oLpewtoop including the rod bundle section. See 


R? 
Appendix A for a derivation of mass flow rate and linear 
power relationships. For a rod bundle section the pressure 
losses through the spacers adds another component to the 


total pressure loss in the loop. The steady state pressure 


balance will then become 


(3-10) 





AP, = APR + AP. = £ AP 


where f. is a multiplying ratio which accounts for spacer 
pressure drop. It 1s introduced in order to make the cal- 
culation of mass flow rates, pressure drops, and linear 
power output easier. The factor f is found to be approxi- 
mately constant for the range of velocities investigated in 


any given lattice configuration. 


Siestcieueang tue relationships for AP. from Eq. (A-3) 


B 
and AP, from Eq. (A-5) into Eq. (3-10) and solving for the 


mass flow rate gives 
gB0-2 sho Sghoe, eqn, 0357 
e ie 
> | i ’ (3211) 


Cmugezi  “sNeuf 
ep s 


w = ( 


where w is the mass flow rate for a unit cell. 


As shown in Eq. (A-3) the thermal buoyancy pressure gain 





ee 
depends on the linear power and the mass flow rate, and the 
friction pressure drop depends on the mass flow rate. 
These dependencies can be seen in Fig. 3.5. In Fig. 3.5 


Eemee(s-li) is the locus of points along the AP,f. line, and 


for any given qi the corresponding w can be taken from the 





graph. The dashed lines represent constant th - te lines. 
Since 
"a Caenic 
gi = 5 (3-12) 
glee 
ag 


for any given core lattice ORAS any one of the three 
weariab les Si w, At can be found, given the other two. 
solutions for a natural circulation operating condition will 
be the point of intersection of the three lines representing 
AP, for a given Ge AP,f., and At. Substituting Eq. (3-11) 
into Eq. (3-12) 


e gB0°2.D.* “Son B85 (2 GSE c, 1.0 
She = Gus => (——) (—) (3-13) 
0.092u>° “NE - 1 
& Ss Nt 


Equation (3-13) is the relationship for the intersection of 
the AP, f. line in Fig. 3.6 with a given At. The quantities © 


w and AP, are taken into account in Eq. (3-13). 


3.6 Comparison of Linear Power Output in Open-hexagonal, 
Triangular, and Sguare Lattices 


As shown in Fig. 1.2, the number of effective rods per 


cell for the three lattices is 


Oe 


1 


2 
<q 


oa ae ae t 
i) ee ! 


Dypreat matural circulation 
operating condition 





a 

oo 

a 

1@) 

ro 

A 

<J 

w (lb/hr) 
oer oe lh iteiratave grapm of solutions for natural 


Piscutauton O©peravang conditions. The intersection of the 
three lines defines an.operating condition. The AP,f line 
is defined by Eq. (A-6). AP,, lines are defined .by Eq. (A-3). 
iiemouisines sare soverned by Eq. (A-11). 
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open-hexagonal peo 2 

triangular = 4 WZ 
u = = 

square n,. uf 


For this comparison it will be assumed for a first order 
approximation that the fluid properties of the coolant, 

At, Ze fos and Ny are equal for the three lattices. The 
three lattices will be compared by taking a ratio of the 
linear power outputs of the three lattices based on equal 
fuel volume fraction and equal pin diameters. The variables 
will be the equivalent diameter, Dos of a flow cell, the 
cross-sectional flow area, S, of a cell, and the number of 
effective fuel rods per cell, ni As shown in Appendix B, 
for equal Vi the hydraulic diameters of the three lattices 


£ 


are equal. By definition 











De 4 xX cross sectional area _ 4S — 
e wetted perimeter P j 
or | Ds 
S = ‘ 
4 
Therefore, 
TAD, 
=) ? 
2 
TAD. 
5S, = A 
8 
TaD 
2 = ae 








Using Eq. 


we obtain 
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(3-13) for q', 
LoS RaE 
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0.092°°*Nne 
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gBo ZD 


Heo oeiG 1.555 Cc 
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q (3-13) 


the following ratios of linear power output: 


1.8) O25 avs) 51D 


(S" “s 
1.8 02555 
rie.) 
1g 


(L/n)y 


ERs) eo, 
Gyn, Jk 


(S 


Ors 515 ie DO 


(1/2) 
OS SS) 5 


((maD,/2) 77° (2)) 
al 


ERS See 


((nav,/8)**8 (1/2) (1/0.5) 


8 (2) 0-9? 13) 335 


ib) 


((naD_/2) 7° (1/2) 


7) 


= 


1.8 Oe 555 ees 5 5 
((7dD./4) ( 


that the linear power output of the three lat- 


tices under natural circulation is the same for equal fuel 


volume fraction and equal pin diameters. 


iiss alo rect 


consequence of the fact that the hydraulic diameters are the 


same. 


Tt is also easily shown that power removal capabilities 


fide tTOnecedmeonveccte1on laminar or turbulent flow conditions 


are also the same for the three lattices (see Appendix D). 


In Chapter 4 some typical values of fluid and lattice 
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parameters will be used to calculate the pressure drops and 
mass flow rates for a range of linear power output. Using 
typical core and primary system parameter values, a q' for 
natural circulation will be calculated to see how this value 
compares to current linear power ratings for forced convection 


& 


cooling. 
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Ghapter 4 


Example Calculations of Mass Flow Rate, Pressure 


EE OE li 


—e SS See eS 


ae ~<Discussion 

In Chapter 3 the relationships for the mass flow rate, 
the pressure drop, and the linear power output per cell 
were derived. In this chapter typical values for lattice 
parameters will be used to solve for the mass flow rate, 
the pressure drop, and linear power output for an open- 
hexagonal lattice and for a triangular lattice of equal fuel 
volume fraction, Vis 
in Chapter 3, the linear power output for triangular, open- 


and equal pin diameter. AS was seen 


hexagonal and square lattice arrangements are equal for 
equal fuel volume fraction, for equal core temperature drop, 
and for equal fluid properties in the three lattices. This 
equality is based on the condition that the pressure drop 
through the flow cells of the three lattices is equal since 
eoolant velocities are equal in the three lattices. This 
condition is not strictly accurate because of the differ- 
ences in the spacer arrangements. These differences intro- 
duce only small second order errors and can be neglected 
Mimasiiesy order calculation of ne relative linear power 
Output. However, since they cannot be neglected in more 
pmeerseucalculacions O: the mass eater rate and linear power 
output of any given lattice arrangement, we will include 


er relations Of | these effects and tales Ooravalues Of mass 
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flow rates, velocities GidmeelesprGosure Groep tcr a range 

of linear power output from one to twenty kilowatts per 

meat for the different lattices. 

4.2 Calculation of Flow Parameters and Theoretical Spacer 
Pressure Drop in the Open-hexagonal Flow Cell 
hemerecoutiec il oa and 1.2.a for definition of the 

parameters required to calculate the expression for the 


Vi of an open~hexagonal lattice. The values of the para- 


meters used to calculate (Vi) are dy = 0.300 in., 6, = 0.040 
mie, and Py, = 0.340 in., from which 
TIJ3 a 
OO), = == ES" = oe: Cay 
9 Py 


The value for dy is about as large a value as is being 
considered for present LMFBR's fueled with uranium dioxide. 


The value for 5, is about as small a pin to pin Spacing as 


h 
is being considered. 

The following values will be used in calculating mass 
flow rates, pressure drop, and linear power output for an 


open-hexagonal lattice flow cell: 


as number of pins per cell = 2, 
ns = cross seevtonal area of a fuel pin = ran /4 
2 


= 0.00049 ft, 
d = assumed spacer web thickness = 0.00125 ft, 
n = number of spacer segments per cell = 9 


1 = length of one spacer segment = p,v3/3 = 0.0164 ft, 


ee a _ 
— 





A 


2 
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ross sectional area of spacers per cell =nl1d 


II 


D Ss s 
SOOO OhOU ft. 
= total cross sectional area of cell = 3 3p6 /2 
= 0.00209 ft 
= flow area away from the gridded GeeaOn. = An ~ nyAy 


0.0011 rt*, 


II 


flow area in the gridded region = Ay - A. 


=0.000918 ft (see Fig. 3.2 for definition of flow 
regions), 


Tc = A,/A, = 0.834. 


Equation (3-11) can now be used to calculate the variation 


pie tlow rate with power output, 


1.2.1.3, a 


2 
Bp ZePe ee a Oliesiey g 


W = (7 a) 9 (Sith) 
0O.,092u Byte | 


with the following additional constants: 


g LB 52 WO? Pape 
B = 1,6 x 107* Sp7t 
0 = 50 lb/ft>, 
n = 0.5 1b/ft-hr, 
c = 0.3 BIU/1b-°F, 
2 2 ) 
D. = (3/3p, - Ta5)/rd, = 0.0281 ft, 
S = A, = @COL1 > 
ie = 2 rods, 
1 See Lo 
3 = 1 (see Eq. (3-10) for a definition of f,). 





NOE 

Thus, for example, for q’ == 10 kw/ft = 34130 BTU/hr-ft, the 
average mass flow rate for an open-hexagonal lattice flow 
cell is ® = 602 lb/hr. Table IV.1 lists values of mass flow 
rate, coolant velocity in the ungridded region of a flow cell 
and corresponding Reynolds number for the open-hexagonal 
lattice. The coolant flow is seen to be turbulent for the 
range of linear power output Conmcidereds 

The following computations are sample Canewrac.ons OF 
the average coolant velocities in and away from the gridded 
region of an open-hexagonal flow cell for the g' = 10 kw/ft 
case. These velocities will be used in PaUeesUolse Wa ee erelaye nels 
Gem cue to the fuel pins and the spacers in the open~hex- 
mena flow cell. The quantity V, is the average coolant 


velocity avay from the gridded region: 
(Vp, = P/p( Ay), = 10,966 ft/ar . 


The quantity V, is the average coolant velccity in the 


eridded regions: 


(V5), = /p(Ag), = 132164 ft/nr 


The quantity Re, is the Reynolds number for the coolant 


away from the gridded region: 
(Re, ),, = Dop(Vy)p/t = 302815 . 


The quantity Re, is the Reynolds number for the coolant 


in the gridded region: 








25,678 







Ef AGH. 


28, 455 


| 29,677 








32,887 
33,840. 
34, 748 
35,614 

| 36, 4d 

| 38, 009 | 

2 


39, 466 


eee 





_o~ 
(Re,), = DP (V5) ,/u is 36,990 


Brom fig, 3.3 for © = 0,834, K, = 0.14 and K, = -0,02. 
Using Eqs. (3-3) and (3-6), the value for the total flow 


cell pressure drop is calculated: 





ie 
(Ap_ - AP_) = AP -+- AP (3-3) 
Jt E RI oar S 
0,09211071 99-81 -8 
(AP, - OP.) = ———p-p— HF) + 
oe 2, Onpe 
p(A,/A5) VS 
2 —=(K + K )n (3-6) 
Oe eC Ces 
Cc 
a :. -6,.1,8 ~6,/2 
(AP, APA a = 1,28 x 10 Vy (O25 ©1052 x 10 Vi 
: (ALE 
(AP, - APL.) EecO.lo/it -. 


: 18 
y ‘1 Se ed 5 
ieamlesnhy.> lists yalues of AP. (PRE h? ((4P Ji ne and 
(AP 


aS AP), for a range of linear power output. Since 


Reynolds numbers and velocities are equal for triangular and 


Gpen=nexagonal lattices, SP,'s will be equal for the two 


R 
lattices, 


4,3 Calculations of Flow Parameters and Theoretical Spacer 
mecescUne cop an tune Triangular Flow Cell 
timer mvOmla Cooma se ang L.2.b for definition of the 


parameters required to calculate the expression for the Ve 
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Table IV.e 


rop in the Open-hexagonal Lattice Flow Cell 











(APg OPA) 
(1b/£t") 





(OP Jenn 


APRf/f ip 
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AP 
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—/4— 
of a triangular lattice. Using tne same fuel volume fraction 


and fuel pin diameter as were used for the open~hexagonal 


mavLLLCe, 
(vi), = (2)? = 0.4 (4-2) 
Vv! = = Oe ae -?2 
ca 293 Dy 


Petvine for the pivtvch and the pin to pin spacing gives: 
on = OME ajo, 
a. = Gi, IWer ayes 
The following values will be used in calculating mass flow 


rates, pressure drop, and linear power output for a trian- 


gular lattice flow cell: 


a ie 

A, = 0.00049 ft’, 

d_ = 0.00125 ft, 

ne =u uae 

0,02 fe; 

A. = 0,0000375 ft*, 

Ap = J3pe/4 = 0,00052 ft°, 

A, = 0.000275 £t°, 

Ay = 0.000237 ft* (see Fig. 3.2 for definition of 


flow regions), 
T = Af A, = 0,865. 
Equation (3-11) can again be used together with the following 


Seq ulowalk constvancs: 


8 


ll 


Hea Gea ad © Pts fae 


em | = 
1.6 x 100° Cr 


& 
B 


3 


I! 


9. 


4 5— 


0 = 50 1b/re?, 
u =N0nS 1b/ft-nr) 
a ORCEEUU Le —, 
(Z./N.) = 0.5s 
2 2 
D, «+ = (2V3pE - Td, )/ad, = 0.0281 ft, 
S = A, = 0.000274 et’, 
al = eat 
f = 4. 
Ss) 


Thus, for example, for g' = 10 kw/ft = 34130 BIU/hr-ft, the 
average mass flow rate for a triangular lattice flow cell 
few 150.5 lb/hr. 

Table IV.3 lists values of mass flow rate, coolant 
velocity away from the gridded region, and the corresponding 
Reynolds number for the triangular lattice. Note that the 
corresponding coolant velocities and Reynolds numbers are 
the same for the open-hexagonal lattice since the fuel 
volume fractions are equal. Note also that the mass flow 
rPaves'for the triangular lattice flow cell are one fourth 
of the corresponding mass flow rates for an open-hexagonal 
lattice flow cell. The differences are due to the differ- 
ences in flow area, Ass per cell and in the number of 
effective fuel rods, n» per cell. The mass flow rates 
per unit flow area, W/D ot? are equal hOmeune uo lLaccvices. 

The following computations are sample calculations of 


the average coolant velocities in and away from the gridded 
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Mable IV.3 


Flow Parameters for the Triangular Lattice 


=e er is ee A EE CE SE | etd 
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region of a triangular flow cell for the q' = 1O kw/ft case. 
diese velocities will be uged in finding the pressure drop 
Omen lo the fuel pins and the spacers in the triangular 
lattice flow cell. The procedure is the same as that used 
in finding the pressure drop in the open-hexagonal lattice 


nllow cell in Section 4.2, 


NG. S66nnt/ hr, 


(Fy), = ela.) 


II 


w/p(A,), = 12,701 ft/nr. 


Note that (V,), = (V,),, since the fuel volume fractions are 


equal, but (V5) 4 is slightly less than (V5), since the grid 


configurations are different for the two lattices, 


II 


(Re, ), = Dp(V,),/A1 = 30,815 


(Reg), = Dp(Vs),/u = 35,688 


muonmeute, 3.3 for UT = 0,865, Ko == (0) SAl2e share! Ka = 30), 08s 
Again using Eqs. (3-3) and (3-6), the value for the total 


flow cell pressure drop is calculated: 


6—2 


(4,095) + 0.036 x 10. Vi» 


a 
(2 = APL). ieee x 10 V, 


30.54 lb/ft? . 


| 


(AP ~ APL), 


S eae De) A 
The lattice pressure drop ( Pp ~ 4P,),, and ( Ee 


are very close in value. This is due to the offsetting 


> Bale 


effects of the values of the f/f,,,, factors and the values 
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of the spacer pressure drop on the values of (4P, - AP.,) 


Ze 

Saee CwOo Lattices. Table 1V.4 lists values of AP 
ft 

(AP_s re ((AP_)e des and (AP, <2 UNE 


cir 
linear power outputs. 


R? 


=) for a range of 


4,4 Empirical Spacer Pressure Drop for the Triangular and 
Open-hexagonal Lattices 


The method of de Stordeur for finding the empirical 
spacer pressure drop in a triangular lattice was presented 
in Section 3.4. Inserting the following values into Eq. 


(3-7), here repeated, and solving for ((AP.) on) ee 


c_V<Sp 
AP =—~2— , (3-7) 
2 2h, 


8 


ll 


h.16 x 10° ft/hr*, 


0,000274 £t°, 


kK 
H 


Le (OSi ae Aone 


ll 


0,0000375 ft*, 
= 50 10/ft?, 


Ogos eee 
! 
>< 
] 


2.33 (see Fig. 3.4 where Re. = d_Vp/u = 1,506), 
gives 

2 
WAR ye Seno ho Gtc-erid |, . 


Using the theoretical spacer pressure drop calculated for 


faveineawbearelatoelce Im Seetion 4.3: 


2D : 
— Aaa ey sikeas a 





~H9O = 


mae Ly, tt 


Pecocire UrOpein the Triangular Lattice Flow Cell 










OE Maree 
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VYherefore, 


Caen 
| (Cia ES = OP Jen oe eH A = 3.56((4P_),), : 


The pressure drop for the open-hexagonal lattice is 


0.09210 1p Oy f 








Oita at eee 2 ) 
e Bo cir 
See) o(A,/8,)°¥ 
S om) A ee c% 1. K,)n, ; 
(AP Jen” 2B, : 


(AP, - AP Jay See 9( 0.25) = 46.3 1y/ft* 


omeculating the factor, is as gives (see Eq. (3-10)): 


ay Ab + (AP_) 
Gan _ R cir s‘em _ 1.96 . 
S aver 


Table IV.5 lists values of V,, Re, Ci, ((AP_) i). 


(OP Dende (COP Jog /(AP, Dende and (aP_).), used in 
finding Ces for a range of linear power output from one 


to twenty kilowatts per foot. (AP and (APf/f are 


Rh cir h 
ROvnGern apie Wyse, The factor (ib Jes is relatively con- 


stant, with an average value of 1.95. The quantity 


(AP,, - AP = (APf oy is also tabulated in Table IV.5. 


ah 
For more accurate pressure drop calculations, Pathe 


should be used in Eq. (A-8) for finding the mass flow rate 
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Weabilkes iV. 5 


Flow Parameters and Pressure Drop in the Spacer Region 
fOrmeniemlmrameular Lattice Flow Cell 
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of the open-hexagonal lattice flow cell. Coolant velocities 
and cell and spacer pressure drop would be recomputed to’ 
find a corrected (foe This process would be continued . 
until succeeding Ces were essentially constant. 

The factor, Pas for the triangular lattice is calcu- 
lated by the same procedure as was used to calculate (fy. 
The average value of (f_) is in ee ine factor, (f_),> is 
not equal to Cia) since the spacer pressure drop in the 
vom lavcices are Ggifferent, a condition due to the differ- 


ence in spacer configuration in the two lattices. 


4.5 Linear Power Output for the Open-nexagonal Lattice 


Using Eq. (3-13), 


2 lige Be: 
Meee) 4s SG mE SS 5.1.0 
eee) (8), (3-13) 
0,092 ° Nf. n, i 


the average linear power for the open-hexagonal lattice 
temecteulaged Usine the following constants: 


SAIS se lee eae 


g 

B a 16 se ee Oe 
9 = 50 1b/rt?, 

UL On lb/tt-hr, 

C, = 0.3 BIU/1b/°3; 


Ll kw = 3413 Bru/hr . 
In terms of several system variables, the linear power 


output can be expressed as 
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aie (ee) 


7 Z Lon eee At pol 
at 22M (ae) (299) | 


S 
4 E 100 nt 


ee eo ; oe : . 
where At is in “F, iD) Pomme tesso 26 in Square inches, 


pam 2s in feet. Inserting the following typical values _ 
into Eq. (4-3): 
gee -TOne cee. the satio of the source/sink 
elevation difference and the effective primary 
circuit length), 


Me =S 4OO ie 


S =e (0), Wea in. (the cross sectional flow area for 
the open-hexagonal flow cell in these calcu- 
lations), 


L. a Lo Gis | 

D), = 0,337 in. (the equivalent diameter for the 
open-hexagonal flow cell in these calculations), 

ns = 2 rodsy 

1 — dy enue, 


the average linear power output for the open-hexagonal 


lattice (and for the triangular and square lattices) is 


qi, = 0.873 kw/ft . 

The above value for g' is low compared to the required full 
power SMnpuineol a pracuical reactor, which averages around 
10 kw/ft. On the other hand, 8.73% of full power is a 
respectable rating in terms of decay heat removal. Sugges- 
tions for increasing the value of the linear power output 


achievable with natural circulation flow are made in 


emaocer oO. 


EY 
The mass flow rates, w, coolant velocities, V, and 


friction pressure drop, AP calculated in the preceding 


R? 
example were based on an ae 1, Table IV.6 lists the mass 
flow rates per unit flow area, W/K 3142 and the total friction 
pressure drop, AP. tos calculated using (f ody = I AS)5 creel 

ie) iin teem sil), The tube bundle friction 


pressure crop, AP ps was calculated from 





APR = i 26 x as yess 3 
where 
_ Te) 
V. = 
1 @ 
pA, 


Figure 4,1 shows the natural circulation operating 
conditions for the open-hexagonal and triangular lattices. 
eae a At = 400 OF, the linear power output of the tyo lat~- 
tices vary by the factor ‘51a aed 


MA cle 


in units of W/A 


Since (A. th = 


—nz 


and w, iw, the horizontal axis of Fig. 4.1 is 


SALI for ease of comparison of the perfor- 
mance of the two lattices. Figure 4.1, which is based on 
Paemvoluece im table 1V.0, shows that for equal total friction 
pressure drop and equal At, the triangular lattice has a 


slightly higher linear power output than does the open- 


hexagonal lattice. 
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Table IV.6 


Mass Flow Rates per Unit Flow Area and Total Friction 
Pressure Drop as Functions of a Calculated t. in the Open- 
hexagonal and Triangular Flow Cells 
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Chapter 5 


Equivalent Diameter for the Open-hexagonal Lattice 


em te eS ee 


5.1 Definition of Equivalent Diameter 

Whenever flow in noncircular ducts is considered, it 
is necessary to employ in the calculation of the Reynolds 
number and friction factor some characteristic dimension. 
In circular ducts the characteristic dimension is the diameter, 
D. In noncircular ducts the concept of equivalent diameter, 
oe (also referred to as hydraulic diameter) applies when 


Pfemauct Joes not appreciably vary from circular. . The 


equivalent diameter is defined as 


sone 1 
p = 4% cross sectional flow area _ 4A (5-1) 


< wetted perimeter Ps 


The equivalent Ate te ig the characteristic length which 
would give the same ratio of fluid passage volume to wetted 
eed tor Circular and noncircular ducts. 
Bpeze COnGitions for a Laminar Flow Region in the Open- 
hexagonal Flow Cell 
In a noncircular duct the concept of equivalent dia- 
meter applies for turbulent or laminar flow. However, the 
values for the Reynolds number and friction factor are 
different for laminar and Phe onl Seer Low. In a unit open- 
Hescsgencaiecel la (terer to Fig. 1.2.a) the flow passage is 


approximately circular; therefore, the concept of equivalent 
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diameter should apply. In Chapter 4 the Reynolds number and 
friction factor for the open-hexagonal cell were calculated 
assuming fully turbulent flow in all parts of the cell. 
These numbers were computed using the average cell rate of 


flow which was turbulent. For turbulent flow 


- -0.2 
fo = 0.046Re ‘ 


However, it is possible that the flow is laminar in some 
sections of the flow cell, namely, in the narrow flow area 
between two adjacent fuel pins. If this is the case, then 
the Reynolds number and friction factor are different in 
the laminar flow area from the Reynolds number and friction 


factor in the rest of the flow area. For laminar flow 


_ -l 
fo = 16Re i 


In order to determine if the flow in the region between 
adjacent fuel pins is fully laminar, a rough calculation of 
the laminar sub-layer thickness around the fuel pins will 
be made. The pin to pin spacing is 0.040 inches; therefore, 
if the laminar sub-layer around each pin is as much as 0.020 
inches, the region may be considered to have laminar flow. 
If this be the case, the open-hexagonal flow cell can be 
divided into two channels, one having turbulent flow and 


one having laminar flow. This division of the flow cell 


assumes no inter-channel mixing and little or no heat trans- 


fer from the laminar region to the turbulent region. This 
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assumption is not good if chee de a significant degree of 
Mieersenwannel Coolant mixing, which in fact probably occurs 
in the present case since the channels are not isolated. 
This assumption breaks down also if the heat transfer from 
channel to channel is good. 

If it is found that open-hexagonal cells can be divided 
into two channels, the coolant flow rates would have to be 
recomputed using two flow regions per cell in order that 
the power output could be recomputed. Since the pressure 


drop through the subassembly, AP will be the same for 


SA’ 
all flow regions, the coolant velocity in the laminar channel 
will be less than the coolant velocity in the turbulent 


channel in accordance with: 


(AP ca) tam ~ (AP oa) mur 
yah ent 
= D.2 earn ~ Af5 2 ) mur 5 


maometals last relation it follows that 


“nena “ature 0 (5-2). 


SiemElemtriction tactor £or E£Uily laminar flow is greater 
Pie eiCmerTCLlon tactor: £or fully turbulent flow in a 
ETmOOLM mipe Or duct, and the laminar region D. Te Ehejoletop<ik— 
mately less than or equal to the turbulent region Doe 
Calculations of the equivalent diameters for the two regions 
in the .cpen-hexagonal flow cell. 


showed (D.) eo) 


- 
Lam e° Tur 
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<2 Eien 9 < Ww where W fom eiemiass 
ge * 1G ; Lam 


Since V 
am Ala sig 


Lam 


flow rate in the laminar channel of the flow cell, and Onur 


is the mass flow rate in the turbulent channel of the flow cell. 
If it 1s desired to maintain the core At at a certain amount, 


+ < rE ley 
then (Wr om Omar?) Woy where Ww, is the flow rate based on 


a Single region of fully turbulent flow in the entire flow 


col, The sum of w and w moles e tiene ts because 
Lam aye sx 


Oy am in the laminar channel is less than Onur? Since oo and 


At are constant and (w,_, + W_,,) < then Sop will have 


am Tur Sie 


to be decreased from the value calculated in a single region 


(turbulent) flow correlation (Gnot) sr = BONE thus degrad- 


ing the power output of the entire core. 


5.3 Calculation of Laminar Sub-layer Thickness in Turbulent 
Flow 


The thickness of the laminar sub-layer around the fuel 
pins will be calculated to a first order approximation using 


a The laminar momentum 


equations given by Schlichting 
sub-layer thickness, Sur for turbulent flow in pipes is 


ouven by Schlichting's Eq. (20.15a) as 


18. 7N) 
64 7, (5-3) 
Piemiateeronme velocity, V,, 1S given by his Eq. (20.8) as 


Vy = (T/o)*/* = (0.03955 (/8y ta Ay 1/2 , (5-4) 


where 


i” is the shear stress at the wall,- 
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is the average velccity in the pipe, 


& | 


Teeuoecnemiinemattc £rScosity of the fluid, 

ad is the pipe diameter. 
Using the following values from the example calculations 
in Chapter 4:3 


qd = (Don Oncol ft 


n= V1 =O oGGe lft /nr 
v = p/p = 0.01 £t?/hr 
egeneriction velocity is v, = 598 ft/hr, and the laminar 


momentum sub-layer thickness 64 ~ 0.001 in. The thermal 
boundary layer thickness on the other hand is Son = ki. 
Using the following approximate values for liquid sodium: 


k 


40 BTU/hr-ft- °F, 


h = 20,000 BrU/hr-£t*-°F, 


the thermal boundary layer is Say =O 002 Sine 

These values of by and Sih are considerably less than 
0.020 inches, the value necessary to establish a laminar 
Flow regicn between adjacent fuel pins. Therefore, the flow 
in all regions of the open-hexagonal flow cell can be 
considered to have turbulent flow. Consequently, the open- 
hexagonal flow cell can be considered to have only one flow 
region, and an equivalent diameter, defined by Eq. (5-1), 
can be calculated for the entire cell. 


The above calculations of Oy and San were only approxi- 


mate. A more accurate computer calculation using a code 


=o 


(13) can ve made 


Ben as General Electric's Velvet-II Code 
of the velocity profiles around the fuel pins. However, 
there 1S no reason to expect that more sophisticated cal- 


culations would alter the above Genecius ions. 
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Ciepter 6 
summary and Conclusions 


6.1 Summary 

Tne work reported in the preceding chapters has compared 
the open-hexagonal, triangular, and square lattice config- 
urations under natural circulation flow with respect to 
leakage reactivity effects in a coolant cell upon sodium 
voiding, mass flow rate and pressure drop in the coolant 
channels, and linear power output. The basis for comparison 
Pimeeees three latvvices was So icaie GOOG (Suze and volumetric 
composition in the cores. The open-hexagonal lattice was 
found to be preferable from the standpoint of leakage reac- 
tivity upon sodium voiding because of its larger coolant 
Channel nuclear hydraulic radius. For the same cell bulk 


temperature rise the mass flow rates per cell flow area, 


Yooii/Acell? 
flow cells, The pressure drop through the lattices were 


were found to be equal for the three lattice 


approximately equal; the small differences in pressure drop 
were due Bostic Cine temeces in "Spacer oat lon! The 
linear power outputs of the three lattices were found to be 
WUeenot thiceednee lattices assuming equal core pressure 
drop. Example calculations of the above quantities were 
presented. It was also found that the fluid flow and heat 
Meeaiouereenmaracteristics Of the Chree lattices were equal 
leer ehoreceagsecirculation for beth laminar and ayyabwsleane 


Miowwiinallivasthe applicability of Che equivalent diameter 
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‘concept for use in calculation heat teansfer and iluid flow 

characteristics of the open-hexagonal flow cell was examined. 
Bievas found Chat the conventional expression for equivalent 


diameter can be used for the open-hexagonal lattice. 


6.2 Discussion 

The open-hexagonal lattice was found to be superior 
to the triangular and square lattices with respect to leakage 
meaculvVivy upon sodium voiding. If more sophisticated 
calculations confirm that the added negative leakage com- 
memeno Or the Sodium void reactivity effect in the open- 
hexagonal lattice is indeed several dollars, there would be 
a definite advantage from a safety standpoint in employing 
this lattice in either natural or forced convection LMFBR's. 

However, from the standpoint or hydraulic and thermal 
performance examined in this study, the open-hexagonal lat- 
ificemvas Net Lound to be Superior to the triangular or 
Square lattices in any way. The triangular lattice would, 
in fact, be preferable to an open-hexagonal lattice since 
more fuel pins can be packed into a given area of core while 
still maintaining prudent pin to pin spacings. Assuming the 
present FFTF pin spacing of 0.0625 inches for a triangular 
lattice and equal pin Memetene of qd, = 0.300 inches for 
both lattices, the fuel pins in the open-hexagonal: lattice 
would have to be touching per Eq. (B-1) in order to retain 
an equal fuel volume fraction in the two lattices. If a pin 
tomo espacing for the open-hexagonal facGkree were to be 


held at 0.040 inches, the Vp as calculated in Appendix C 





GIS) 
would have to be held to around 0.306 without considering 
stainless steel swelling and to 0.288 with consideration 
Given to stainless steel swelling. Larger pin sizes for 


the open-hexagonal lattice would allow larger V but the 


f? 
increase in Vp would not be significant if the pin to pin 
spacing was held at a reasonable value (6, < 0.020 inches). 
This can be seen from the decreasing slope cf the Ve versus 
a Curve Of fig >l. 1. Current and future design values for 
Ve are found eee Oe 6. 1, 

Since the natural convection performance of alJj. three 
lattices is equal in the hydraulic and thermal aspects 


examined, improvements in the linear power output achieved 


by all three lattices can be suggested using Eq. (A-10), 


eee se lO 

oi gBep Z_D > nn oO ye 1 nO 

ol = es oa) ie 52D (Aas . (A-10) 
OSE et ae ny. 1 


In Chapter 4 a value for Z/N the ratio of the source/ 
sink elevation difference and the effective primary circuit 
Hengzth, of 0.5 was employed. This ratio can be increased. 
Using values from the FFTF design proposals for a 400 Mut 


uvrpr''4), 5 typical value for N_ is 


e 
AP ii 
i ee, ay an lp Sie 
e Core : 
AP 
core 
where 
AP = HUG jets! 


pump 
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AP ope = 100 psi 


Leore 7~ 8 tt 


By employing a Z, of up to say 24 Feet, the ratio ib AE 
could be made as large as 2.0. For a value of Ze Ne = 20), 
the linear power output calculated in Chapter 4 would be 
increased by a factor of (MOONS) 2 = eee eine Chapter 
4 the value of the core temperature rise was 400 °F, If 
this value could be raised to 450 °F, the linear power out- 
put would be increased by a factor of (450/400) 1+??? =e 
A core At of 400 °F is about as large a value as is prac- 
tical. A larger value of At may produce a thermal shock 
upon reactor scram which the primary system could not with- 
stand. The linear power output could also be increased by 
Spacing the fuel pins farther apart while increasing pin 
diameter. The wider spacing would increase the equivalent 
diameter, D,, and coolant flow area, S, thus increasing q' 


Ripe) 2>) Z Th aaa & 


Since q' = £((D The core 
temperature rise, At, would also increase with this change 
in spacing through the relation qn, = a een At: However, 
the spacing contribution probably would not increase q' by 
much. In Chapter 4 the core length, 1, was four feet. The 
tength 1 could possibly be reduced by as much as two feet, 
thus raising q' by the factor (4/ayr© = 2, Combining the 
above suggested contributions for increasing q' gives a 


total factor of about 5.2. With these contributions q'! 


Would equal 4.5/| kw/ft instead of 0.873 kw/ft as calculated 
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in Chapter 4. However, +.54 kw/ft is in reality a peak 
merue Lor tie natural circulation case. It is less than 
50% of the desired 10 kw/ft average and less than 30% On 
the 16 kw/ft peak eommene Aoetwmanadsradi@al peaking factors 
of about 1.25 and 1.26) needed to match projected 1000 Mwe | 
mereed convection plant operating some aons ae! Unless a 
better combination of the variables discussed above can be 
achieved, full power operation under natural circulation is 
mee practical. 

A summary of the parameters assumed and calculated in 


this study is presented in Table VI.1. 


6.3 jphessiunere Work 

Since the leakage reactivity effect in the open-hexagonal 
lattice provides a significant additional negative reactivity 
contribution compared to the triangular and square lattices, 
a more sophisticated calculation by Monte Carlo or other 
methods should be made to determine a precise value for the 
sodium void negative reactivity component in the open- 
hexagonal lattice. The model used in this study yielded 
only approximate answers under the assumption of an cnvetath be 
core. An experiment measuring the relative and absolute 
negative reactivity of the lattices should also be consid- 
ered. However, since the reactivity change is on the order 
of one percent of the effective multiplication factor, the 
necessary precision for measuring the reactivity change 


would be difficult to design into the experiment. 
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Models of the three’ lattices could be constructed to 
test the pressure drop performance of the different config- 


Mracbtons. Freon or water could be used as the test coolant. 


6.4 Conclusion 

This study has shown that from hydraulic and thermal 
aspects, the triangular, open-hexagonal and square lattices 
ees ecoentially equivalent under forced and natural circu- 
lation. For natural circulation operation it was shown 
that the linear power output of the three lattices was far 
less than the desired value for full power operation, even 
under an optimum combination of constraints. Therefore, 
memewOk une lattices are practical for full power operation 
under natural circulation. Furthermore, unless the open- 
hexagonal lattice is found to contribute significantly to the 
negative component of the sodium void reactivity effect, 
wiemenianegulear lattice is preferable for forced convection 


operation. 
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Appendix A 
Maes OWerenetatronsnip tor Natural Circulation 
Hg, 2 Reactor 
If z is the elevation above a datum surface (see Fig. 
A.1) at which the pressure is constant, the difference in 
pressure between elevations z. and z, in a static fluid 


Z i 
(no friction nor acceleration) is 


Ey 
( edz _ p&(2Z, - 2p) 


ee eee A-1 
Oo I g Eo ( ) 





Integrating completely around a closed circuit which nas 


the fluid in two legs at different densities gives 


aD iia a 
” zl f P22 + J p,a2), 
fo Ey 
bs 
& J 
AP, = E (P, - P.)dz . 
& 


Introducing the temperature coefficient of volumetric 


EDans Lon, © ¢ 


Ph 7 Po * PB, (ty : te) 


AP, eee se {, 3 Ges (A-2) 





Uisabeves ve] =we (t = t) where q is the heat produced per 


h 
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Cooler 


i Se | (Heat 
Exchanger ) 





Primary Coolant 


(16) 


ioe liernaleemecilataon driving heads. 





~(3- 
miieuecell, and @ is the mass flow rate per cell and 


= 2 Bo SO aa aes Z)/2 


gB.p .qZ 
ee, (A-3) 
Bot yy 


From the detinition of the Fanning friction factor, 


AS Els oan |B, 
i ae ; (A-4) 
ON pu 


ean expression for the friction pressure drop, APR» is 


obtained as follows: 


= 
tT 
2N pw fpf. 


AP. = ; (A-5) 
R gD.p°S” 


where ives Hicwucaevorececouncing Lor Spacer pressure drop, 


is defined by Eq. (3-10), and u=w/pS. Equating AP, and 


AP, and solving for w gives 


o— as 
x gBp-qZ_.D_s 
eee /3 (A-6) 


eN.cpt Bt s 


In this derivation Ne is the total effective length of the 


primary coolant system: 


AP 

N = Ee 

e AP cCOLre 
core 
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where L, is the length of the core. The quantity N, is 


ore 
a measure of the equivalent length of the complete primary 


circuit (core, pipes, and heat exchanger) expressed in 


terms of an extended core of height N and channel diameter 


D,. Substituting the relationship f, = 0.046Re7°"* 


= 0.0461°°*/(p.W/S)°** for turbulent flow into Eq. (A-6) 


gives 
. (A-7) 


Using the relationships gq = nail, where n, is the number 


Semcitective fuel rods per unit cell, and q FS LENG. where 


ae = th ~ t.: 





q' = : (A-&) 


Substituting Eq. (A-7) into Eq. (A-8) yields 


_ e6p°z_pi-?st-Pn gi 0. 357 c At 
qt = ( ms y0-357(-P) . (A-9) 
0.092 Necpts ic 


Solving for q', Eq. (A-9) becomes 


Beer llatrc 


2 
gBp 4D S 


0.092u NAF, n,. 


CLO 
MN 


“_ Mr y0.555(01y1-555(°p)1-0 (a0) 
L 
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Refer to Fig. 3.6, which is reprocuced below as Fig. A.2, 
honmewne paysical illustration of theese Gecmebva. stele. 5 
snows a comparison of the driving forces in the forced and | 
hatural convection cases. 

For natural convection the following observations are 
made: If it is desired to increase q', this is done by 
allowing an increase in se, This in effect increases the 
driving head allowing a greater flow rate, wover Eq. (A-8) 
at which the AP = APR: The combination of an increase in 
both At and w increases q according to the relation 
eect = t.)» and thereby gq! increases for fixed n,, and 
1, since q! = q/n,1. 

Note if At is fixed, i. ©., a maximum allowable core 
temperature drop is seeded - the cell operating condition 
for constant cell parameters n. and 1 is fixed at the inter- 
R? and AP... ice ce r= 
section fixes ™ and hence g and q' through q = DPedio nie tines 


section of the three lines At, AP 


and q! = q/n,1 = We At/n 1 from Eq. (A-10). 





AP 





a. Intersection points given by Eq. (A-}. 





b. Eq. (A-10) yields the gq! which corresponds to the 
Mmiemsecellonmpoimecmonetme three lines; i, e., for con- 


Garaomls Ot AP, = AP, Kacey fixed. 


IHN sa\erear Natural circulation operating conditions. 


-([- 


AP 
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a. Relationship between driving pressure and mass 
iew rave in a forced eonvection system where, for a fixed 
At, q' = we At/n 1, and w depends on the pump pressure. 


AP 





w 


b. Relationships between driving pressure, mass flow 
rate, and core temperature rise in a natural convection 


system where, for a fixed At, q' = Wo At/n1 is found by Ea. 
(A-10), and w is aetermined as the intersections of lines 
oi Ge APB» and AP 2° 


winch COUPartsomeot Graying forces in forced and 
Het wae cCOnvect LON. 





Appendix B 


Derivation of Hydraulic Diameters as Functions of 


[aelietemensnaettine undies cemmmintoandtin tani) = == Sa eee —— a 


A EI a ee ee 


Refer to Fig l,l for finding fuel volume fractions, Ve? 


2 2 
Sida / rd dq 
ee Te) (B-1) 
“7 65p8 8B 38re 9 py, 
T .2 2 : 
a /4 Ta Toa 
2 x one Ne 
Cn = ~t— = a (—-)~ , (B-2) 
sa p58 2)3p, 6 Py 
2 
mam /4 Tt d 
2 
(Vi). = —B— = -(--)*. (B-3) 
q p hy - 


In the derivations of the linear power output capa- 
bilities of the lattices, the assumptions made were that 
the same fuel pin diameter was to be used in each lattice, 
and the fuel volume fractions were to be equal. Hence, we 


obtain the following relations between lattice pitches: 


(Vi)y = (Vi), 5) 
(By? = 2)? | 


33 Py 243 D, 
and 


ps /pe SOB ye (B-4) 





2 SiS 
In like manner 


2 & 
P sq’ Ph = sayz (B-5) 


4 x cross sectional flow area 
mov cetinition De = ae 
wetted perimeter 





b(56py (p,/2) 13 - ends /t) 3 Gve - Ta. 


So, SOE eet 
Pedn = ord 1d ii 
Pp p 
W504 (P,/2)03 - srae/i) 2 pt - 
On) a ee -  (B-7) 
2p p 


Using Eq. (B-4), 


2{3(Sp5) - wae 
(Dh he & Ene = (Dae (B-3) 
p 
U(p2 - ra@/4) Up ~ raé 
(D.) Sq 7 Thy) sq 7 9 ; (B-9) 
wee Td TD 


Using Eq. (B-5), 


4 (393p2) 


(D..) = (Oy (B-10) 


sq” 1d 
p 


The above results show that for equal pin diameters and 
fuel volume fractions, the equivalent diameters of the 


three lattices are equal. 
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Appendix © 
Lattice Parameters 


Table C.1l lists some of the lattice parameters whicn 
are being used in current large LMFBR's and which are 
Prone sed ror iuvure large LMFBR!'s. Table C.2 and Fig. C.1 
show volume fraction of fuel and trianguiar lattice pin 
pitch and gap spacing based on a gap spacing of 0.O40 inches 
for the open-hexagonal lattice (the minimum assumed prac- 
tical) and equal pin diameters for both lattices. 

The equations used in these calculations are Eq. (B-1) 
and Eq. (B-2). The calculation procedure was to select d., 


erdpeusimg bo. = 0.040 inches, compute DP, and (V;), from 


h 


Eq. (B-1). Then taking (Vin = (Vi) and 6, are computed 


Cost f 
from Eq. (B-2) for the same assumed pin diameter. Fig. C.1 
shows that as q., decreases, O, approaches oF with an accom- 


panying drop in Vi. For reasonable pin sizes as shown, and 


EOMSUbaInes Of constant oh and V! 


P between lattices, the 


following oeanne one are drawn: 
Tine MeLOn pm eZap fOr a triangular lattice 
( 0.120 in.) is considerably larger than that for an 
open-hexagonal lattice CO SOMO shine 
Pe iiceniel volumes traction for the lattices with 
reasonable pin diameters (0.250 to 0.300 inches) are 
about 0.45 to 0.47. Significant increases above these 
values require large (and uneconomical) increases in 


pinsdiame ten. 
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Table C.,2 


Mimucmodmi stacrne and pin pitch for triangular 
ana open-nexagonal lattices which have equal 
Pedmetamevers and equal volume fractions of fuel. 
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Seeetucesuel yolume fraction, Ves used in the above | 
calculations was based on a single mueio cell in a fuel 
subassembly and was equal to 1 = Vy (sodium volume 
iene Mem— ating into account the volume or the core 
taken up by the clad, the gap between clad and fuel, 
the subassembly housing, the space between subassem- 
blies, and the control ageciinilies., Shite actual fuel 
volume fraction, Vp» for the whole core is 0.306, as 
will be shown in detail below. 
Ceo eevomuro nC. fOr Carculavion cf the actual core 


Mmicievolume fraction, V The fuel subassembly contains 


ee 
olf fuel pins of 0.300 inches outside diameter and has a 
Vi of Cee (2eeebaccewonmea Unig lenecvh of fuel assembly, the 
inside volume of the fuel subassembly is 217V,, 24 patie ie 
= Se in.°, Dividing the subassembly into six equilateral 
triangles, the height, h, of a triangle becomes one half 

the internal distance across the flats. Using a wall thick- 
ness of 0.150 in. and one half the subassembly to subassem- 
bly spacing of 0.025 in., the subassembly pitch divided 

Sec vOulout eee O,f7>5 in, One equilateral triangle of 
the subassembly has a volume of (32.2) =e. 5 in.”, The 
volume of one triangle equals sph = S(S3)n = v3, = Sao, 
in.” Ooi eee melon! e= 3.04 + 0.175 = 3.215 in. 

The volume each subassembly occupies, Von in the core based 
on the distance from subassembly center across the flats 


to the midpoint betveen subassemblies is Van = 6(Sb'h') 
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Section a. 


Pic Waco ce cuDassemb1y configuration 


2 ; 
= 6 (£(248nt yn) = 2 fBni* =o nils in.?. Using a pin diameter 
BTS) 
peo OOM ma clad thickness of 0.015 in., and a fuel to 
clad radial gap of 0.002 in., the volume of fuel in a sub- 
Be _ POs sO0 = 0.034) ae: 
assembly is aT = ee) OO Fiat 
Assuming there is one control assembly for every ten fuel 
subassemblies, the actual fuel volume fraction in the core, 


e is 


V 12706 
V . = oot = (~——_______—-) = 0.306 
Von 4- Tosa Boon tes 351 5 


This number is low compared to the fuel volume frac~ 
tions cP = 29,2 to 46.7) of the reactors in Tabie C.1. 
Our value of Ve does not take into consideration the 
irradiation induced swelling of the stainless steel sub- 
assembly walls. Because of this swelling the subassembly 
to subassembly spacing must be 0.25 inches iiovueea O10, 05 
iMmenecwu thnerenores i) = 3,04 + 0.275 = 3.315 in., and 
Von = 38.07 in., Consequently, our value of V, is lower: 
va = 0.288. The fuel volume fractions of the pores in 
Table C.1 do not reflect the stainless steel swelling. The 
methods used to calculate the volumetric compositions of 
the cores for the reactors in Table C.1 were not explained 
an the references cited. Therefore, it is difficult to make 


a good comparison of the volumetric compositions of the 


Ae ciaecorece in ta0le GC. and our theoretical core. However, 


if the V, of 0.306 is tco low for practical use, an increased 


ie 
pin diameter is feasible with the use of uranium carbide 
fuel which has a higher melting point than does uranium 
dioxide fuel. Some designs are proposing pin diameters 


of 0.400 inches for future 1000 Mwe LMFBR's fueled with 
uch 154) , 
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- Appendix D 


Se ee ee em ee ee ee 


The average linear power output of a lattice flow cell 


is given by 


usc At 
Cl a ae (oe) 
nt 
where 
w= AV 


Assuming the pumps are moving the coolant at the same total 
rate through the three lattices, and in view of the fact 
that the volumetric compositions of the three lattices are 
equal, the coolant velocity will be the same in the three 
meee S elas Shown in Appendix B that the equivalent 
diameters of the three lattices are equal. The lattices 
have the following numbers of effective rods per flow cell: 
a) aon oar == 2 ove = 1. 


Since 


Di x (wetted perimeter) 


A= F (D-3) 
m 
IDL eeupts) D_1d 
AL = peleray) Pomp (D-4) 
dy B) . 
D.(3r4,,) ; 


oC 
A, a a 


D 1d 
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ge = hy (D-5) 





-9Je- 


iL 
=e ae = 5A, . (D-6) 


Using the above expressions for flow areas in Eq. (D-2), 
which is substituted into Eq. (D-1), gives the following 
relative values for the linear power output of the three 


lattices: 


ay, - pA, Ve At/el - A, /2 a 
— wer... pl. ae ,L 
f ae — 

qt pA Ve At/5l GAL/5 
Gq pA,Ve,At/21 A, /2 , 
qq pA Vc, At/1 A, /2 

a! 

oe i 

sq 


The above ratios show that the linear power output for the 


three lattices are equal under forced circulation for 


f 
The pressure losses through the different lattices are 


equal Vi and or 


the same since 


hipvé 
EF 5) (D-7 ) 





and all the parameters in Eq. (D-7), including f, and D,, are 


B 
Pir euoimttiemcaree Lauthees, The corresponding coolant 





ToIleocneluesemyseare Siown vo De equal in the lattices in 


Mavpleswivy.il and £¥.3. The forced convection heat transfer 


Sectimuetenus are also the same for the three lattices since 


Nu 


n 
hD./k = C, + C,(Pe) ; 


Nu = C, + Co(D.Vpe,/k)” (D-8) 


and all the values in Eq. (D-8) are again equal for the three 
laeci.ces, including the D's as discussed in Appendix B. 

From the above ratios and correlations it is seen that 
the imdcanspOwer OUCOUtL, pressure drop, and heat transfer 
characteristics are the same for the three lattice config- 


ire tLOns. 
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English 


Appendix F 


Nomenclature 
(unless defined as they appear) 


Meaning 
cross sectional flow area 
heat capacity 
Civaite vic 
Pistetnom — accor 
gravitational constant 
heat transfer coefficient 
thermal conductivity 
iUmeroLLCapiom Taccor 
length 
a quantity 
pipe length 


Nusselt number 


average power output cer cell 
average linear power output per pin 
pitch | 

perimeter 


Peclet number 
DigeSouice drop 
radius 

Reynolds number 


cross sectional flow area 
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ides 
rt", in. 
BTU/1Lb-CF 
fee Mls 
et/hr- 
BDU/hpafe = ° 8 
BIU/nr~-ft- oP 
if Were ae 
tay 
kw 
kw/ft 
aide. 
cgi 


1b/ft*, psi 


alitelee 


rt , in. 
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Meaning Units 
O 


axial temperature difference per cell i 
(Ginmic vaso. OUG LES ) 


average coolant velocity ft/hr 
fuel volume fraction =1 - oe ---- 

where Va is the sodium volume 

fraction 

Actual fuel volume fraction in a -—— = 


real core 
height ct 


elevation difference iG 


fraction of delayed neutrons a 


temperature coefficient of Qed 
volumetric expansion 

Divneue sol spacing al 
boundary layer thickness siete 
dynamic viscosity lb/ft-hr 
kinematic a SiGhOOERGey et-/br 
density lb/ft? 
ratio of areas, Ao/Ay re 
average mass flow rate per flow Llb/nr 


cell 


due to buoyancy forces 
equivalent 
empirical 


entrance, exit 
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SUDScripts ' Meaning 
gz grids 
h open-hexagonal lattice 
hy hydraulic 
p pin 
a rods 
R que to friction 
S spacers 
SA subassembly 
' sq square 
c triangular Lattice 
th TieOre vical: 
Al total 
W wetted 
1 away from gridded region 


2 in gridded region 
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